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The critical currents in Josephson junctions of conventional superconductors with macroscopic 
defects are calculated for different defect critical current densities as a function of the magnetic 
field. We also study the evolution of the different modes with the defect position, at zero external 
field. We study the stability of the solutions and derive simple arguments, that could help the 
defect characterization. In most cases a reentrant behavior is seen, where both a maximum and a 
minimum current exist. 
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I. INTRODUCTION 

The interaction of localized magnetic flux (fluxons) 
with defects (natural or artificial) or impurities in su- 
perconductors or junctions has an important effect in 
the properties of bulk superconductors or the behavior of 
Josephson junctions correspondingly Q. The flux trap- 
ping from defects which is of major importance in Joseph- 
son junctions ^J can modify the properties of poly crys- 
talline materials with physical dislocations, for example 
grain boundary junctions \^. In this category one can 
also consider grain boundary junctions in YBa2Cu3 07 
M where the tunneling current is a strongly varying 
function along the boundary. This strong inhomogeneity 
makes them good candidates for SQUID type structures 
[g| . Phenomenologically the current-voltage characteris- 
tics of grain boundary junctions are well described Iq] 
by the resistively shunted junction model . The grain 
boundary lines often tend to curve, while the junction 
is very inhomogeneous and contains nonsuperconducting 
impurities and facets of different length scales ^^ . The 
linear increase of the critical current with length in grain 
boundary junctions with high- Tc superconductors, which 
is a different behavior from the saturation in the inline 
geometry of a perfect junction, can be explained by the 
presence of impurities |l^ . Therefore it is interesting to 
study flux trapping in impurities especially when it can 
be controlled. Modern fabrication techniques can with 
relative ease engineer any defect configuration in an ex- 
tremely controlled way. 

In bulk materials are several types of defects that can 
influence the critical current in high temperature su- 
perconductors like YBa2Cu302: materials. They include 
3d inclusions, 2d grain boundaries and twin boundaries, 
and point defects like dopants substitutions, oxygen va- 
cancies m. For example the homogeneous precipitation 
of flne Y2BaCu05 non-superconducting particles in the 
melt processing of YBa2Cu302; leads to high Jc values 
due to the particle pinning centers pT| . Similar be- 
havior is observed in NdBa2Cu302; bulk crystals with 



Nd4Ba2Cu20io particles [Of. Of interest is also the case 
of the peak effect in twin-free yi23 with oxygen defi- 
ciency. In this case, one sees a linear increase (peak ef- 
fect) of the critical current at small magnetic fields, when 
growth is under oxygen reduction |T3]. For the fully ox- 
idized crystal one expects a decrease. The peak effect 
is attributed to flux trapping. Information on the defect 
density and activation energies can also be obtained from 
the I-V characteristics, as was the case for several types of 
defects which were also compared to Au^ irradiated sam- 
ples with artiflcial columnar defects ||lj| . These columnar 
defects also act to trap flux lines in an YBCO film which 
is considered as a network of intergrain Josephson junc- 
tions modulated by the defects. In this case assuming a 
distribution of contact lengths one finds a plateau in the 
critical current density vs. the logarithm of the field p5[ . 

The study of long size of impurities is going to give 
information beyond theories which concern small ampli- 
tude of inhomogeneities [Q. Also it is possible for a 
direct comparison of the numerical results with exper- 
iments in long junctions obtained with electron beam 
lithography Q. This is a powerful technique which al- 
lows the preparation and control of arrays of pinning cen- 
ters. Another method is the ionic irradiation which pro- 
duces a particular kind of disordered arrays, consisting of 
nanosized columnar defects ||lj,|l3| . The variation of the 
critical current density can also occur due to temperature 
gradients ||19| . 

The activity in the area of high critical current den- 
sities in the presence of a magnetic field is hampered 
by defects due to the difficulty of having a high quality 
junction with a very thin intermediate layer. Thus sig- 
nificant activity has been devoted, since for example the 
energy resolution of the SQUID [E0| and the maximum 
operating frequency of the single flux quantum logic cir- 
cuit pi[ , to name a few applications, depend inversely 
and directly respectively on the plasma frequency Wp, 

1/2 

with LOp r^ Jc ■ The fundamental response frequency 
of Josephson devices, the Josephson frequency loj, also 
depends on the critical current density. On the other 



hand, a drawback is that high critical current densities 
lead to large subgap leakage currents |22| ] and junction 
characteristics degrade rapidly with increasing Jc- 

Variations in the critical current density also influence 
the I — V characteristics introducing steps under the in- 
fluence of both a static bias current and the irradiation 
with microwaves p3| . In that case the variation is quite 
smooth (of sech type ), so that the fluxon and its mo- 
tion could be described by a small number of collective 
coordinates. Interesting behavior is also seen in both the 
static and dynamic properties for the case of a spatially 
modulated Jc with the existence of " supersoliton" exci- 
tations p^,p5| and the case of columnar defects [ p6|j27| 
or disordered defects | [lO| |. 

The trapping of fluxons can be seen in the Imax{H) 
curves where we also expect important hysteresis phe- 
nomena when scanning the external magnetic field. The 
hysteresis can be due to two reasons: (i) One is due to the 
non-monotonic relation between flux and external mag- 
netic field [ p8[ arising from the induced internal currents, 
and (ii) from the trapping or detrapping of fluxons by de- 
fects. The effect of a defect on a fluxon and the strength 
of the depinning field depends strongly in the size of the 
defect, the type of defect and the position of the defect. 
Here we will consider case where the widths of the de- 
fects is of the order of the Josephson penetration depth. 
In this range we expect the strongest coupling between 
fiuxons and defects. We will also consider the case of a 
few defects in the low magnetic field region where pinning 
and coercive effects are important. 

The organization of the paper is as follows. In Sec. 2 
the sine-Gordon model for a Josephson junction is pre- 
sented. In Sec. 3 we present the results of the critical 
current Imax versus the magnetic field of a junction with 
an asymmetrically positioned defect. The variation of the 
Imax and the flux content Nf with the defect critical cur- 
rent density and the position are presented in sections 4 
and 5 respectively. The effect of multiple pinning centers 
is examined in sections 6 and 7. In Sec. 8 we examine 
a defect with a smooth variation of the critical current 
density. In the last section we summarize our results. 



II. THE JUNCTION GEOMETRY 

The electrodynamics of a long Josephson junction is 
characterized from the phase difference ^(a;) of the order 
parameter in the two superconducting regions. The spa- 
tial variation of (f>{x) induces a local magnetic field given 
by the expression 



nix) 



d(j)[x) 
dx 



(1) 



in units of Ho, = 



2TTdX, 



where $o is the quantum of 



flux, d is the magnetic thickness and Xj is the Josephson 
penetration depth. The magnetic thickness is given by 
d — 2Xl + t where Al is the London penetration depth 



in the two superconductors and t is the oxide layer thick- 
ness. The Xj is also taken as the unit of length. The 
current transport across the junction is taken to be along 
the z direction. We describe a 1-D junction with width 
w (normalized to Xj) in the y direction, small compared 
to unity. The normalized length in the a;— direction is 
£. The superconducting phase difference 4>{x) across the 
defected junction is then the solution of the sine-Gordon 
equation 



d^4){x) 
dx'^ 



Jc(a;) sin[0(x)], 



with the inline boundary condition 



dx 



I 



x=±i 



H, 



(2) 



(3) 



where / and H are thejrormalized bias current and ex- 
ternal magnetic field. Jc{x) is the local critical current 
density which is Jc = 1 in the homogeneous part of the 
junction and Jc — jd in the defect. Thus the spatially 
varying critical current density is normalized to its value 
in the undefected part of the junction Jq and the Xj used 
above is given by 



A,/ = 



$n 



27r/io(iJo ' 



where /zq is the free space magnetic permeability. One 
can also define a spatially dependent Josephson penetra- 
tion depth by introducing Jdx) instead of Jg. This is 
a more useful quantity in the case of weak distributed 
defects. 

In the case of overlap boundary conditions Eqs. (0) 
and (0) are modified as 



d^(t){x) 
dx"^ 



Jc{x)s\-a[(l){x)] - /, 



and 



dx 



H. 



(4) 



(5) 



We can classify the different solutions obtained from 
Eq. (0) with their magnetic flux content 



Nf 



1 



^(^a-'/'l), 



(6) 



in units of <I>Oj where (t>R(L) is the value of 4>{x) at the 
right(left) edge of the junction. Knowing the magnetic 
flux one can also obtain the magnetization from 



M^^Nf 



H. (7) 

For the perfect junction, a quantity of interest is the crit- 
ical magnetic field for flux penetration from the edges, 
denoted by Hci- For a long junction it is equal to 2 



while for a short it depends on the junction length. Due 
to the existence of the defect this value can be modified 
since we have the possibility of trapping at the defects. 
For a short junction we have penetration of the external 
field in the junction length, so that the magnetization ap- 
proaches zero. For a long junction it is a non-monotonic 
function of the external field H . 

To check the stability we consider small perturbations 
u{x, t) = u(a;)e'** on the static solution 0(x), and linearize 
the time-dependent sine-Gordon equation to obtain: 



cPv{x) 



Jc{x) cos (j){x)v{x) ~ \v(x)^ 



(8) 



under the boundary conditions 
dv(x) 



dx 



0, 



where A = — i 



It is seen that if the eigenvalue equa- 



tion has a negative eigenvalue the static solution 4>{x) 
is unstable. There is considerable eigenvalue crossing so 
that we must monitor several low eigenvalues. This is 
especially true near the onset of instabilities. 



III. ASYMMETRIC DEFECT 

In the following we will consider the variation of the 
maximum critical current as a function of the magnetic 
field for several defect structures. We start with a long 
(L > Xj) junction of normalized length £ = 10 with a 
defect of length d = 2 which is placed D = lA from 
the right edge. Thus the defect is of the order of Xj. 
We plot in Fig. 2a the maximum critical current Imax 
variation with the magnetic field. The different curves 
correspond to phase distributions for which we have a 
maximum current at a given value of the magnetic field 
H. The overlapping curves called modes have different 
fiux content as seen in Fig. 2b where we plot the magnetic 
fiux in units of <f>o for zero current versus the external 
field. The magnetic flux is only a weak function of the 
external current. 

For the perfect junction there is no overlap in the mag- 
netic fiux between the different modes. In fact each mode 
has fiux content between n$o and (n -I- l)$o and there- 
fore is labelled the (n, n + 1) mode [^. Here in the case 
of the defect the range (at zero current) of fiux for each 
mode can be quite different and the labelling is with a sin- 
gle index n = 0, 1, 2, ... corresponding in several cases to 
the (0,1), (1,2), (2,3),... modes of the perfect junction. 
There are in several cases several modes with similar flux. 
To distinguish them we add a letter following the index 
n. 

The maximum Imax is obtained for mode 1 and the 
increase comes from the trapping of flux by the defect. 
We have to note that the (d, e) part of this mode is a 
continuation of the (a, b) part of mode 0. In both cases 
we have entrance of flux from the no defect part of the 



junction and the instability in the critical current occurs 
when (p{—£/2) — tt. Here and in the following we will 
take this to mean equal to tt modulo 27r. For the maxi- 
mum current (at H < 0) the equation is iJ — 1/2 = —2. 
This can be understood from the pendulum phase di- 
agram, where the (j)x = —2 is the extremum slope, and 
thus the relation Imax — 4-I-27J holds. For the (&, c) part 
of mode the flux enters from the right where the defect 
is. This reduces the critical current compared to the per- 
fect junction mode [g8|,|9|. Note that 0-mode has its 
critical current Imax peak slightly to the left of iJ = in 
the Imax vs H diagram, and to the left of Nf = in an 
Imax vs Nf diagram (see Fig. 11a). Also in the absence 
of current, reversing the direction of H only changes the 
sign of the slope d(j)/dx, but the phase difference (in ab- 
solute value) at the two ends will be the same. Thus the 
mode at / = extends for —1.6 < H < 1.6. This is 
not clearly seen due to curve overlapping in the left side. 
Comparing the Imax for the modes 1 and —1 we see that 
Imax{^) > Irnaxi—^)- In both cascs & fluxou (or anti- 
fiuxon) is trapped in the defect. The major difference in 
the Imax comes mainly from the phase distribution which 
in the mode 1 case leads to a large positive net current in 
the undetected side, while in the —1 mode the net current 
in the undetected side is very small. 

For the mode 1, at iif w and zero current the insta- 
bility happens due to the competition of the slope of the 
phase at the defect center and at the right edge, while 
at the other end a,t H = 2, the field at the defect center 
becomes equal to the external field applied at the bound- 
aries and there is no such competition. In this case the 
instability sets in due to the critical value of the phase 
at the undetected boundary (i.e. (l)x{—i/2) — 2). The 
situation is analogous for the mode —1. For iJ « the 
instability sets in due to the depinning of the antifiuxon 
while for H = —2 due to the critical value of the phase 
at the free defect part of the junction. For the mode 
we have no fiuxon trapping at the defect, even though 
the instability at the two extremes with H — ±1.6 at 
zero current is caused from the tendency to trap a fiuxon 
or antifiuxon correspondingly at the defect. At higher 
values of the magnetic field (|7?| > 1.6) we have stabil- 
ity for a range of non- vanishing current values as will be 
discussed below. Thus this value can be considered as 
the minimum value for the introduction of fiuxons in the 
junction. Let as remark that for the perfect junction, 
or a junction with a centered defect, the corresponding 
values for fiuxon introduction would be equal to 2. Thus 
there is a decrease of the critical field as the defect moves 
away from the center. For the mode a centered defect 
would have no influence on the solution. 

The results for the maximum current are in agreement 
with the stability analysis. In Fig. 2c we present the low- 
est eigenvalue Ai for the different modes in zero external 
current / = as a function of the magnetic field H. The 
sudden change in slope for the modes —1, 1 is because at 
that point a new eigenvalue becomes lower. The Ai is 
positive denoting stability and becomes zero at the criti- 



cal value of the magnetic field, where a mode terminates. 
The symmetry about zero magnetic field is due to the 
symmetric boundary conditions for / = 0. Change of 
the sign of H changes the sign of the phase distribution, 
but the cos <j) in (pi) remains unchanged. This symmetry 
is being lost when a finite current is also applied. Also 
there are solutions (not presented in the figure) for which 
the stability analysis gives negative eigenvalues i.e insta- 
bility. These solutions may be stabilized when we insert 
multiple impurities. 

In Fig. 3a we specifically draw only the 1 mode, to 
be discussed in more detail. Here we changed the proce- 
dure, in searching for the maximum current. Up to now 
we followed the standard experimental procedure, i.e. we 
scan the magnetic field and for each value of H we in- 
crease the current /, starting from / = 0, until we reach 
the maximum current. Here we consider the possibility 
that for / > there is also a lower bound in the value of 
the current for some values of the magnetic field. This 
requires a search where we vary both H and / simulta- 
neously. Thus we see that for H < there is a lower 
bound given approximately by the line H + 1/2 « Hd , 
where Hd « is the critical value of iJ at / = 0, for 
which we have depinning of the trapped fluxon. Over 
this curve the slope (j)x at the right end (near defect) is 
kept constant and equal to Hd and above this line the 
fluxon remains pinned and it should be stable. This line 
ends at H = —1, since in that case the extremum value 
(j}x = —2 is reached in the left end. Increasing now in 
that range of H the bias current we find that also the 
Imax curve extends further to the left. The equation for 
this line is approximately given by if — 1/2 = —2, with 
an extremum at 4>x{—i/2) = —2. Thus the instability on 
this line arises from the left side (far from the defect). It 
extends up to 77 = —1 for a long junction and joins the 
other Hue H + 1/2 = Hd- 

The above calculations were done for a long junction 
so that the fields at the two ends do not interfere. For 
shorter length however the two ends feel each other and 
in that case the two instabilities are not independent. 
This means that the tail of the defect free side field will 
compete with the slope of the trapped field. Then the 
two lines H + 1/2 = Hd, and H — 1/2 = —2, end before 
they meet (at H « —1) at a cutoff magnetic field. Also 
for short junctions we expect the straight lines to have 
some curvature. A similar discussion holds for the right 
end of the 1 mode. Again there is a lower current (pos- 
itive) bound given by if — 1/2 = 2 due to instability at 
the left end {(f>x{—2) = 2), and an upper bound given by 
H + 1/2 = Her, where Her ~ 2.8, due to fluxon depin- 
ning. On the same diagram, we show the lower bound 
for negative currents. Thus we see that there is strong 
asymmetry for positive and negative currents. Remark 
that for negative currents the mode 1 is very similar to 



the mode (1,2) with no defect |28|| 

right boundary is determined by 

undetected side. The left boundary is again very close 

because Hd w 0. So an interesting effect of the defect 



. This is because the 
an instability at the 



is that we have this strong asymmetry for positive and 
negative currents. 

We would get a similar picture if we considered the —1 
mode. In fact we get the same curves (as for mode 1) if 
we put / — > — / and H => ^H . This is consistent with 
the —1 mode shown in Fig. 2a. The discussion can also 
be extended to the other modes. In Fig. 3b we show the 
result of a similar scan for the mode, but for the sake of 
shortness we will not discuss the —1, —2, 2 modes. In any 
case when the number of fluxons increases one must rely 
on numerical calculations rather than simple arguments. 



IV. VARIATION WITH THE DEFECT CRITICAL 
CURRENT 

In the previous section we considered the case of a 
microresistance defect. With present day masking tech- 
niques we can also consider any finite critical current 
(lower or higher) in the defect. This situation can also 
arise very often in junctions with high critical current 
densities, where small variations in the thickness can 
create strong critical current density variations. Thus 
for the previous asymmetric defect configurations we will 
study the effect of the defect critical current density in 
the magnetic interference pattern Imax [H) ■ We will con- 
centrate on the and 1 modes. 

(i) mode 0; 

In Fig. 4 we see the Imax{H) variation for the mode 
for decreasing values of the critical current density from 
id = 2, to j^ = 0. Let us discuss first the case for jd <'^- 
For the perfect junction where jd = 1 we have a sym- 
metric distribution about H = 0. As we decrease jd the 
flux content of this mode (and the extremum H) is sym- 
metrically reduced (see Fig. 2b). It is not apparent from 
the drawing, due to the superposition of several curves 
on the left side of the diagram, but as expected the range 
of the magnetic field is symmetric about ii = at zero 
current. The corresponding Iraax{H) curves, however, 
are not symmetric. The right hand side of the curves is 
displaced towards smaller critical fields with decreasing 
jd- This means that the critical field at i = to in- 
troduce a fluxon from the ends is decreased due to the 
existence of the defect which acts with an attractive force 
on the fiuxon. The curves are linear and can be approx- 
imated by the equation I{H) = 4 — 2{H + 6Hc), where 
SHc is the decrease in the critical field Hco and depends 
on jd- A similar decrease happens for negative magnetic 
fields where the defect tries to pin an antifluxon. Even 
for higher currents the right side critical flcld is deter- 
mined by the tendency of the defect to attract a fluxon. 
The left hand side, however remains rigid (but is shifted 
along the line). This is due to the entrance of magnetic 
fiux from that part of the junction where there is no de- 
fect. The instability in the critical current occurs when 
{(j}{—£/2) = tt) for every value oi jd- From the pendulum 
phase diagram which is the classical analog of the Joseph- 
son junction, the extremum occurs at 9a;^(— 2) — ~2, or 



H — 1/2 = — 2 which is the equation for this triangular 
side. At near zero current the critical field is influenced 
from the attractive action of the defect. At low cur- 
rents and extreme negative magnetic field the Imax curve 
shows a re-entrance behavior so that it is not stable at 
low and high currents, but only for a finite intermediate 
range of current values. This way we reconcile the dif- 
ferent origins of the instability mechanisms 0(— //2) = tt 
at high current and the defect influence discussed for the 
right hand side of the mode. 

For j(i > 1 we see an increase in the Imax, while the 
critical magnetic field at / = remains almost constant 
at about Her ~ 1-9. The instability at that point is due 
to the trapping of flux in the region between the positive 
defect and the right edge of the junction. The field for 
that is expected to be near _ff = 2 if the right undefected 
part is of length of the order of Josephson length. Thus 
it is the same value for flux penetration from the perfect 
junction edges. It will vary weakly with jj^. 

(ii) mode 1: 

The mode 1 in the perfect junction has a full fluxon for 
magnetic field H = 0.07. The phase distribution is about 
^ -K where the energy has a minimum. At the end of 
this mode at H = 2.07 where two fiuxons have entered 
the junction the phase changes from (f){—£/2) = — tt to 
(f>{i/2) = Svr. When the defect is inserted this mode 
is significantly modified due to the fiux trapping in the 
defect. 

In Fig. 5 we see the magnetic interference pattern for 
this mode for different values of the defect critical cur- 
rent density. For < jd < 0.7 the Imax vs H curves 
are displaced downwards, and a fluxon is trapped in the 
defect. We notice that all the curves for jd < 0.7 have 
the same critical magnetic field iJ = 2 for / = 0. This is 
because at this end of the mode, at / = the instabihty 
arises at the side with no defect where the phase reaches 
the critical value <j) = tt (modulo 2n). Of course as dis- 
cussed in the previous section we have a reentrant behav- 
ior above H — 2. At the other end for small magnetic 
field the instability is due to depinning of the trapped 
fiuxon. For 0.7 < jd < 1.0 the defect can trap the flux 
only for H < Had, where the value of the Had depends 
on the defect critical current jd and in Fig. 5 it is shown 
for jd = 0.9. Notice that for this value of jd the fluxon 
is very weakly trapped, and the untrapping process hap- 
pens slowly over a range of magnetic field values. For 
H > Hcd the fiuxon has moved away from the defect, 
and for this weak defect the junction does not feel it. 
The critical current goes abruptly close to the curve for 
the perfect junction. We conclude that the behavior of 
the junction for values of j^ close to jd = 1 is determined 
by the ability of the defect to trap one fluxon. This can 
be seen also from the change in the lowest eigenvalue vari- 
ation with the external field H, at values of the critical 
density jd > 0.7, in Fig. 6. 

For jd > 1 (thin lines in Fig. 5) it has a similar form 
as for jd = 1, i.e. there is no fiuxon trapping. Again, as 
in the mode, the Her at / = stays around 2.0 and is 



again due to the trapping of flux in the right edge. 

In Fig. 7 we present the evolution of Imax with the 
defect critical current density jd for a magnetic field H = 
1.5. For this value of the magnetic field there are no 
solutions with trapped fiuxons for jd > 0.83. The lowest 
eigenvalue at / = becomes zero at this point. For 
jd > 0.83 and H > 1.5 there are solutions which are 
not trapped. For these solutions the maximum current 
coincides with the one of the perfect junction and there 
is a discontinuity in the curves. Notice the point at jd = 
1.0. In the same figure we also show the magnetic fiux at 
/ = and at Imax which is almost constant as a function 
of jd as expected, with small difference between the two 
different current curves. 



V. VARIATION WITH THE DEFECT POSITION 

In Fig. 8a we see the evolution of the critical current 
at zero magnetic field as we move the defect from the 
right edge of the junction D — to the left edge where 
D — 8. The position is measured from the edge of the 
junction to the nearest edge of the defect. We examine 
the several modes separately: 

(i) mode 

For this mode and for / = 0, we are able to find solu- 
tions for all the defect positions. As we can see in Fig. 8b 
the corresponding magnetic flux at Imax is slowly chang- 
ing and equal to zero when the defect is in the junction 
center. But when the defect is placed close to the ends 
the magnetic fiux at the maximum current deviates from 
zero. The critical current for this mode is symmetric for 
defect positions about the junction center, and has its 
maximum value when the defect is at the center. This 
is because at that position it does not infiuence the so- 
lution at the edges which is very close to the undefected 
case, while near the center the phase is almost zero. But 
when the defect comes close to the junction ends the de- 
fect cuts into the area by which the current fiows, and 
the critical current is reduced. 

For even smaller distances D = 0.2, and D — there 
is a jump to solutions which correspond to a current, 
which is much higher than that of the mode for nearby 
D values. This is because the defect cuts negative current 
regions and for this position we have an increase of the 
critical current. In fact these solutions (see -I--I- symbols 
in Fig. 8a) are very close to the solutions of a perfect 
junction within the undefected area, except that now the 
defect at the edge can give contribution to the flux but 
no contribution to the current. Thus the flux is much 
higher than that of the mode and it approaches that of 
mode 1. Nevertheless these points should be considered 
as a separate mode. In fact they are part of a branch 
(crosses). In these distances there are no other modes 
for H — 0. Similar results were obtained by Chow et al. 
|pO| where they attributed this enhancement in the Imax 
for small distances to a self field which was generated by 



the current, penetrating into the defect and resisted any 
further penetration of field. To overcome this resistance 
it was necessary to apply a higher current. But they do 
not distinguish between modes with different flux con- 
tent, and their evolution with the defect position. 

(ii) modes I, —1 

For these modes we do not have solutions for all the de- 
fect positions at / = and H = 0, but only in the range 
1.4 < D < 6.6 as seen in Fig. 8c where the lowest eigen- 
value is plotted as a function of the defect position for 
the different modes. The curves for the 1 and -1 modes 
coincide, while the mode shows a change of slope corre- 
sponding to the last two points { D — and 0.2 discussed 
above) which belong to another curve. Mode —1 has a 
trapped antifluxon in the defect. When the defect is to 
the left (4.0 < D < 6.5), then the instability in the cur- 
rent of mode —1 at i/ = is created at the right-end 
of the junction when the phase reaches the critical value 
(j){l/2) — TT. This instability occurs for currents which are 
less than those necessary to unpin the antifluxon. Notice 
that increasing the current there is no competition with 
the slope of the antifluxon trapped in the left end. Thus 
at this point (for 4.0 < D < 6.5) the maximum current 
is very close to the undetected junction mode 0, except 
that in this case Nf w — 1 is close to an antifluxon. At 
the other end {D < 4.0) the instability for mode —1 is 
caused by the depinning action of the applied current, 
which takes now much smaller values (close to zero) be- 
cause of competition with the pinned fluxon. The phase 
distribution at the defect free end is that expected for 
ff — and / close to zero. The mode 1 with a fluxon 
trapped has a symmetrically reflected (about the center) 
form in Imax vs D and the instability for D > 4.0 occurs 
at the left end of the junction, which is the opposite case 
of mode —1. The eigenvalue becomes zero at the posi- 
tions D — lA and D — 6.6. The Xi{D) curve coincides 
for the modes 1 and —1 due to the fact that the phase 
distributions for the same D for these modes are sym- 
metric about X = L/2, and the cos(j){x) that enters the 
eigenvalue equation is the same. 

In the rest we examine the variation of the critical value 
at which the instability sets in, as we scan the magnetic 
field in the positive (negative) direction Hcr{Hci) for zero 
current, for the different modes, as a function of defect 
position. This instability can be attributed to the pin- 
ning, or the depinning field or to the critical value of -p 
at the defect free edge, depending on the particular mode 
that we are considering. Explicitly for the mode the in- 
stability in the Hci{Hcr) is due to the pinning of a fiuxon 
(antifiuxon), respectively. In this mode the defect has no 
infiuencc for positions close to the center as seen in Fig. 
9a and 9b. However as we move the defect close to the 
edges the pinning field Hci[Hcr) is reduced in absolute 
value because it is easier to trap a fluxon (antifluxon). 
For the mode 1 the Her is constant for all defect posi- 
tions. This is due to the fact that at J = it is the phase 
distribution at the undetected edge of the junction that 
determines the instability. Notice that due to the reen- 



trant character the critical magnetic field takes higher 
values at larger bias currents which vary with defect po- 
sition. The Hci curve depends on the phase distribution 
near the defect and therefore is strongly defect position 
dependent. For the mode —1 the picture is reversed com- 
pared with the 1 mode. In this case the Hd is constant 
while the Her varies with position. Note that in this 
mode the depinning of an antifluxon is the reason that 
causes the instability at Her- 



VI. TWO SYMMETRIC PINNING CENTERS 

As noted defects (with jd < 1) or inhomogeneities in 
the junction can play the role of pinning centers for a 
fluxon. In this section we discuss more precisely the ef- 
fect of multiple pinning centers on the magnetic inter- 
ference patterns Imax{H) and the flux distribution. The 
pinning effect of the Josephson junction has also been 
analyzed in [pl^,p2| , by using a simple mechanical analog. 
The analogies of the mixed state of type II superconduc- 
tors and vortex state of the Josephson junction has been 
discussed in these references. In Fig 10a we present, as 
an example, the critical current Imax versus the magnetic 
field for a junction which contains two defects of length 
d = 2 placed symmetrically at a distance D = 2 from 
the junction's edges. We examine the following modes 
grouped according to flux content: 

(i) modes 0, Oa 

These modes have magnetic flux antisymmetrical 
around zero field, as seen form Fig. 10b where the mag- 
netic flux is plotted versus the magnetic field. At / = 
and magnetic field H = —0.7, the Oa mode contains one 
fluxon trapped in the left defect, while an antifluxon ex- 
ists at the other part of the junction. As H increases 
towards 0.7 the picture changes slowly, so that the anti- 
fiuxon is pinned in the right defect. The stability analysis 
shows that this mode is unstable. We remark that there 
are also other unstable modes near zero flux, which we 
will not present here. For example there is another un- 
stable mode with the same flux as Oa but a much higher 
critical current (the same as the mode). Mode has 
phase distributions which are similar to the correspond- 
ing mode of the homogeneous junction since it has no 
trapped flux in each defect. 

(ii) modes 11, Ir 

These modes have magnetic flux close to unity, and 
are both stable. For the mode Ir one fluxon has been 
trapped to the right defect, and in the mode 1/, the vor- 
tex is trapped in the left defect. Due to the symmetry 
this mode has the same magnetization as the mode Ir, 
but the critical current is reduced. The phase distribu- 
tion for the modes Ir, and 1^, at zero current are related 
by (/)i;(a;) = 2tt — (J3ir{—x). The maximum field H — 1.9 
(at 7 = 0) for both modes is determined by an instabil- 
ity at the defect free side. At the other extreme there 
is a competition at the fluxon side between the applied 



field and the field created by the pinned fiuxon. Thus the 
critical field a.t H = —0.62 can be considered as a coer- 
cive field and below this value the fiuxon gets unpinned. 
The two modes have characteristically different currents 
and this depends on the current through the fiuxon free 
defect, since the pinned fiuxon itself gives no major con- 
tribution. Thus the maximum current is much larger for 
the Ir mode. The opposite would be true if we look for 
negative currents. There are also the symmetrically sit- 
uated modes that correspond to an antifiuxon in the left 
or right defect, which are not shown in Fig. 10a. The 
respective fiux is antisymmetric with H around H ~ 0. 

In Fig. 10a we also show the mode 2 with fiux around 
2 fiuxons. Several unstable modes are not shown, for 
the sake of clarity. Their analysis however, can show 
the connection between different modes, while a defect 
in the correct place with proper characteristics can sta- 
bilize these solutions. We conclude that depending on 
the positions where the vortex is trapped we may have 
modes with the same magnetic flux content, but differ- 
ent critical currents. Also due to soliton localization on 
the defects, we may have stable states with magnetic fiux 
close to unity, for zero magnetic field. These states to- 
gether with the one existing in the homogeneous junction 
form a collection of stable states in a large H interval. We 
must comment here that states with unit fiux, for zero 
magnetic field {H — 0) exist in the homogeneous junc- 
tion, as a continuation of the stable (1) mode to negative 
magnetic fields, but as we found in a previous work p8[ , 
are unstable. So we may argue here that the presence of 
defects stabilizes these states. 

In comparing the results for one (Fig. 2a) and two 
defects (Fig. 10a) we see some similarities and differ- 
ences. In the case of two defects new modes appear but 
also the region of stability of the equivalent modes is dif- 
ferent. This is more clearly seen in Fig. 11 where we 
plot the Imax vs Nf for both cases. This presentation is 
useful since the Nf is a nonlinear function of H. This 
plot (Fig. 11a) is a combination of Figs. 2a and 2b. We 
should point out that the maximum peak in the current 
in both case comes due to the trapping of a fiuxon in the 
defect at the right side. The maximum of 0-mode is very 
close in both cases and this happens because this mode 
does not involve fiuxon trapping. The Ir mode for the 
two defect case is very close to the 1 mode of the single 
defect, since in both cases there is a fiuxon trapped in the 
same side. In the two defect case we see an enlargement 
of the region of stability so that the modes overlap. The 
thin continuation lines in modes and 1 for the single 
defect are in the reentrant region of fiux as discussed in 
section 2. 



VII. SYMMETRIC DISTRIBUTION OF PINNING 
CENTERS 

In this section we study as an example the case where 
a junction of length £ ^ 14.2 contains three defects of 



length d = 2, and the distance between them is 2. The 
length was augmented, so that we keep the same width of 
the defects when we increase the number of the defects, 
since we saw that the width of the order d = 2, gives 
the possibility of fiuxon trapping and increased maxi- 
mum current when the defect is situated asymmetrically. 
We will study the phase distribution at / = and try 
to extract information about the critical field values and 
magnetization. We find the following modes grouped ac- 
cording to fiux content: 

i) modes 0, 0/, Or, Oc 

In Fig. 12a we present the critical current versus the 
magnetic field for the modes with magnetic fiux around 
zero (see Fig. 12d). This is indicated by the sym- 
bol. There are four modes belonging in this category, 
which are stable. The solutions for the mode are sim- 
ilar to the homogeneous junction mode 0, with no fiux 
trapping in the defects. The only difference is that the in- 
stability in the critical field occurs when the phase at one 
edge, reaches a value, which is smaller (due to pinning) 
than the corresponding value for the undetected junction, 
which is (j){—i/2) < tt. The same was true for the two 
defect case. Mode Oc has the maximum critical current 
Imax = 5.08 for H = 0. One antifiuxon is trapped to 
the leftmost defect, one fiuxon to the rightmost, and the 
phase in the center defect is constant. The trapping at 
the edge defects leads to a positive current distribution 
between them, for this particular length, and enlarges the 
maximum current. The same type of mode was not found 
for the two defect case (with a shorter junction length), 
and we conclude that the extra defect along with the in- 
creased junction length stabilizes this solution. For the 
mode 01 one fiuxon is trapped in the left defect where the 
phase changes about the value (jj — t:. The antifiuxon is 
distributed at the other two defects, where the phase is 
about the values 3tt/2 (or tt/2), and we have a cancella- 
tion of the positive and negative current density in this 
region. Similar for the mode Or the fiuxon is trapped 
to the right defect, and the current is distributed with 
opposite sign to the other two defects. These modes are 
similar to the Qa mode for the two defect case. 

a) modes 11, Ic, Ir 

In Fig. 12b we see the maximum current versus the 
magnetic field for the modes with magnetic fiux around 
Nf = 1 (see Fig. 12e). There are three modes with 
fiux close to Nf = 1 each of which corresponds to the 
trapping of one fiuxon in one defect. In the mode Ic the 
fiuxon is trapped in the center defect. In the mode 11 
(Ir) it has been trapped in left (right) defect. Due to 
the symmetry the lowest eigenvalue, and the magnetic 
fiux coincides for these two modes, but as we showed in 
the previous section, their critical currents are different, 
depending on the tunneling current distribution in the 
region with no trapping. The Ir mode corresponds to a 
higher critical current. 

in) modes 2, 2a, 2b 

In Fig. 12c we see the maximum current versus the 
magnetic field for the modes with magnetic fiux around 



Nf = 2 (see Fig. 12f). Only the mode 2 corresponds 
to stable solutions. There we have two fluxons trapped 
in the side defects. In mode 2a one fluxon is trapped in 
the right defect, while in mode 26 this trapping occurs in 
the center defect. We conclude that distributed pinning 
centers are more effective in trapping the vortex, and lead 
to an increased critical current. Some conclusions will 
continue to be valid for larger number of defects where 
we keep the defect width and separation fixed. In that 
case we also expect the results to change significantly 
when there is either a periodic array of defects, where we 
expect higher fluxon modes to give the highest current 
peak pi. 



VIII. DEFECT WITH A SMOOTH VARIATION 
OF CURRENT DENSITY 

Up to now we considered defects with abrupt changes 
in the local critical current density and the question arises 
whether the abruptness of Jc variation is crucial in the 
significant change in Imax for the n — 1 mode. We will 
see that similar effects exist for smoother variation, where 
again the fluxon pinning is an important feature. For this 
reason we chose a single defect at the junction center with 
a smoothly varying critical current density given by 



Jc{x) = tanh 



-{X - Xq) 



(9) 



where the defect is centered at Xq, and the width is de- 
termined by /i. In Fig. 13 we show the results for the 
case xo — 7.6 and /i = 2, which can be compared with 
the results of the asymmetric defect in Fig. 2a. For the 
modes shown the curves are very similar and thus we see 
that the main results survive since the defect strengths 
are similar. Of course there is a quantitative difference. 
But most of the stability criteria described earlier are still 
valid. 

In Fig. 14 we consider the effect of the form of current 
input and compare the case of inline with overlap for 
a smooth defect situated at the center of the junction 
i.e. Xo = with fi = 0.5. In Fig. 14a we present the 
Imax for inline boundary conditions, and we show only 
the —1,0,1 modes. The mode is not influenced at all 
from the defect since all the phase variation is at the 
boundaries. There is a strong similarity with Imax for 
the 1, and —1 modes. The reason is that in these cases 
there is a trapped fluxon or antifluxon at the center and 
at zero current and magnetic field the phase variation dies 
out at the boundaries. Thus when increasing the current 
at H = towards Imax we have the same situation at the 
boundaries as for the mode and the instability happens 
at close Imax values. Of course due to the pinning, the 
fluxon content is very different from the mode. The 
— 1 and 1 modes have an enhanced Imax and the small 
difference in Imax from the mode is attributed to the 
small inffuence of the trapped ffuxon to the boundaries. 



Let us remark that a similar situation was seen in Fig. 
8a for the square well defect, when the defect position is 
at the center for H = {). By comparing with Fig. 9 the 
Hci and Her values we see close agreement with the case 
of id = in the defect. These results could change for a 
smaller length junction or if we move the defect towards 
the edges (as seen in Fig. 8a). 

For the same defect we also investigated the effect 
of the overlap current input, where the current is dis- 
tributed along the whole junction. In Fig. 14(b) we 
present the maximum current per unit junction length 
versus the magnetic field, and it should be compared to 
the inline case in Fig. 14(a). We see a significant change 
for the —1 and 1 modes. Of course at / = both current 
inputs give the same solutions, but Imax is much smaller 
for the overlap boundary conditions. This is from the 
fact that due to the applied current the fluxon is pushed 
against the pinning barrier until it is overcome at the crit- 
ical current. In the absence of applied current the phase 
at the defect center is ^(0) = tt, while the application of 
the current pushes the fluxon to the edge of the defect 
which is taken to be near the point where the curvature 
of the defect critical current distribution changes sign. 
So we can consider in this case this maximum current as 
a measure of the pinning force. 

In Fig. 15a we plot the magnetic flux N f at zero cur- 
rent versus the magnetic fleld H for the inline case. The 
lowest eigenvalues for the different modes versus the mag- 
netic field are seen in Fig. 15b. For a homogeneous 
junction the mode is the only stable state available at 
H — 0. However in the problem we consider here, the 
mode 1(— 1) exists and it is stable for H = and corre- 
sponds to the localization of the soliton (antisoliton) in 
the inhomogeneity. For these modes we have pinned flux 
at ff = 0, with (j){0) = TT, and ^ = 2. In Fig. 16 we 
show the evolution of and ■^, for the mode 1, as we 
change the magnetic field at / = 0. Near H = —1.9 the 
fluxon content is near zero and for H < —1.9 an insta- 
bility sets in due to the depinning of the fluxon. This is 
because the slope at the pinned fluxon competes with the 
opposite slope tried to be imposed by the external nega- 
tive magnetic field at the boundaries. At the other end 
the flux is equal to 2, and the instability sets in when 
at the boundaries approaches n (or odd multiplies). The 
range of H values for the 1 mode, when the defect is at 
the center is significantly broadened and gives a corre- 
sponding range for the flux of two fluxons. Usually each 
mode has about one extra fluxon and in particular for 
the perfect junction it contains only one extra fluxon. 
This is because the defect is at the center and far enough 
from the edges where the magnetic field is applied and 
therefore even for negative fields there is no significant 
competition, with the field at the defect center. This is 
especially true when the distance of the defect from the 
edges is grater than 2Aj. When however the defect is 
near the edge the instability sets in before we cross to 
negative magnetic fields. 



The maximum current Imax for the mode is greater 
than in the modes 1,-1, but is reduced compared with 
the Imax for the mode in the homogeneous junction, in 
zero field. In [p3[ they approximated this reduction in an 
analytical calculation using a delta function for the defect 
potential, and they found Almax = — /i/2L « 0.02. In 
p3[ they arrived at the analytical result, by minimizing 
the fluxon free energy, for the maximum overlap current 
versus the magnetic field H, for these modes, which is 
a good approximation of the numerical solution we con- 
sider here in the limit L ^ 1. 



IX. CONCLUSIONS 

In several applications it is desirable to work in an ex- 
tremum of the current for a region of the magnetic field. 
This can be achieved by the appropriate distribution of 
defects so that the negative lobes of the current distri- 
bution in the junction due to the fluxons are trapped in 
the defect with no contribution to the current. Of course 
if the defect is isolated (far from other defects or the 
edges) we expect zero contribution to the current. Due 
to the effect of the applied current and magnetic field at 
the boundaries in certain cases we can obtain positive 
current lobes outside the defect. In several cases in sec- 
tion 3 this was the reason for the increased current. Be- 
cause the control of magnetic field is very easy compared 
to other system parameters (like temperature, disorder, 
etc.) the measurement of the effect of magnetic field on 
junction behavior, provides a convenient probe for the 
junction. The calculation of the Imax can characterize 
the quality of the junction or verify the assumed dis- 
tribution of defects when they are artificially produced. 
The spatial variation of the critical current density on low 
Tc-layered junctions, and high Tc grain boundary junc- 
tions can be directly imaged with a spatial resolution of 
Ifim using low temperature scanning electron microscopy 
{LTSEM) [ pipq] . Information on smaller scale inhomo- 
geneities has to rely on the magnetic field dependence of 
the maximum tunneling current Imax- 

The purpose of this paper is the consideration of large 
defects in order to study the interaction between fiux- 
ons and defects and give estimates of the coercive field 
for pinning or depinning of a fluxon from a defect. The 
region of consideration puts us far from the region of 
perturbation calculations and is amenable to direct ex- 
perimental verification since it is easy to design a junc- 
tion with the above characteristics. The defects infiuence 
strongly the low fluxon modes. At high magnetic flelds 
larger than the depinning field of a single fluxon we ex- 
pect only minor effect and fluxon trapping. Of course 
for a large number of defects interesting behavior can be 
obtained. [Q^ The interaction between fluxons in the 
few defect case also assists to overcome coercive fields 
and untrap fluxons. The results of two trapped fluxons 
in the two defect case show that the fluxons are strongly 



coupled and one cannot consider an exponential interac- 
tion type potential between the fluxons. Also the critical 
current in a long junction, cannot be calculated as the 
Fourier transform of the spatial distribution of the critical 
current density Jc{x), at least for weak magnetic fields. 
For strong magnetic fields, where we have the field pen- 
etrating uniformly the junction, as is the case for short 
junctions, we recover the diffraction like pattern. 

In summary we saw that the bounds of the different 
modes determined by the stability analysis depend on 
two factors: (i) the instability at the boundaries away 
from the defect when (l)^ reaches its extremal values equal 
to ±2, and (ii) the instability due to the pinning or depin- 
ning of a fluxon by the defect. If the junction is near one 
end then we saw that both criteria play a role in deter- 
mining the instability, independently in different areas. 
In general, however, there will be coupling between de- 
fects and the edges (surface defects) especially in the case 
of multiple defects. Defects also introduce hysteresis phe- 
nomena which are weaker in the case of smooth defects. 
We also saw that due to fluxon trapping in general we see 
a reentrant behavior, i.e. there are regions of magnetic 
field for which there is both an upper and a lower bound 
on the maximum current. We also find that due to the 
pinning of magnetic flux from the defect there exist ad- 
ditional stable states in a large interval of the magnetic 
field. The abrupt change in the critical current density is 
not crucial for the trapping. Similar results are expected 
from smooth defects, with quantitative differences. The 
above results can be checked experimentally since it is 
easy to design a junction with a particular defect struc- 
ture, using masking techniques. In fact a few parameters 
or characteristics could give at least partial information 
on defect properties. In particular the measurement of 
Her or Hci can give some information of the defects near 
the edges. Also one can imagine the situation where we 
scan locally with an electron beam affecting thus the lo- 
cal critical current and observe the variation of the Imax 
as we increase the heating. Once a fluxon is trapped we 
can decrease the heating (or increase jd) and observe the 
variation of Imax ■ Thus one can have pieces of informa- 
tion to put together in guessing the defect structure that 
might fit the whole Imax pattern. The extension to many 
defects requires considerable numerical work. It is hoped, 
however, that some of the stability criteria will still be 
useful. 
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FIG. 1. The geometry of the junction. The dark shaded 

region marks the defect in the intermediate layer, i is the 

junction length and D the separation between the left edges 

of the defect and the junction. 
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FIG. 3. (a) Critical values of the bias current as a function 
of the magnetic field. The solid line is drawn for the mode 
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FIG. 10. (a) Critical current Imax and (b) magnetic flux 
Nf, versus the magnetic field H, for the different modes, for 
a junction of length £ = 10, which contains two symmetric 
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FIG. 11. Critical current 7c, versus the magnetic flux Nf, 
at the maximum current for the different modes, for a junction 
of length i — 10, (a) for the asymmetric defect case, and (b) 
for the two symmetric pinning centers of length d — 2. 
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and (c) 2, 2a, 26. for a junction of length t = 14.2, which 
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The corresponding magnetic flux is presented in (d), (e), (f). 
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FIG. 13. Inline critical current Imax versus the magnetic 
field i/, for the modes 0, 1 and 2, for the junction with an 
asymmetric defect but smooth variation of the critical current 
density.xo = 7.6 and \i = 2. 
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FIG. 14. Critical current Ima,^ versus the magnetic field H^ 
for the different modes, for (a) inline current and (b) overlap 
current, for the junction with a centered defect, and smooth 
variation of the critical current density. 
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FIG. 15. (a) Magnetic flux at zero current as a function 
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in Fig. 14. (b) The corresponding evolution of the lowest 
eigenvalue Ai for the different modes. At the end of each 
mode the Ai vanishes. 
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FIG. 16. The evolution of (bix) and 



J^^' , for the mode 1, 



as we change the magnetic field at / = 0, for the same type 
inhomogeneity as in Fig. 14. Numbers are H values. 
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